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Abstract 

We investigate the simulation of fermionic systems on a quantum computer. 
We show in detail how quantum computers avoid the dynamical sign problem 
present in classical simulations of these systems, therefore reducing a problem 
believed to be of exponential complexity into one of polynomial complexity. 
The key to our demonstration is the spin-particle connection (or generalized 
Jordan- Wigner transformation) that allows exact algebraic invertible map- 
pings of operators with different statistical properties. We give an explicit 
implementation of a simple problem using a quantum computer based on 

standard qubits. 
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I. INTRODUCTION 



Because of recent exciting algorithms, like the factoring algorithm of Shor [J1J and the 
search algorithm of Grover [jj, that solve difficult problems on a quantum computer using 
algorithms that would be impractical on a classical computer, it is easy to overlook that the 
original proposals for quantum computers were for the purpose of solving quantum physics 
problems ||. People like Feynman || focused on the extent to which such a computer 
could imitate a specific physical process, suggesting in part that quantum problems were 
inherently too complex for a classical computer 0]. 

The obvious difficulty with deterministically solving a quantum many-body problem 
(of fermions or bosons) on a classical computer is the exponentially large basis set needed 
(i.e., the dimension of its Hilbert space grows exponentially with the number of degrees 
of freedom). Exact diagonalization approaches (e.g., the Lanczos method) suffer from this 
exponential "catastrophe". Viewed the other way around, this basis set scaling is what 
restricts today's classical computer to simulating only small quantum computers. This 
point seems indisputable, but should not be taken as proof that quantum systems cannot 
be simulated on a classical computer. By the same token, the recent claim |4| that quantum 
computers can simulate all quantum systems efficiently lacks explicit and detailed algorithms 
for specific problems, and lacks a generic model of quantum computation including the 
unitary maps (quantum gates) that can be physically implementable. Even if a quantum 
computer existed, some interesting quantum problems, such as finding the ground state of a 
general quantum Hamiltonian, do not yet have efficient quantum algorithms. Finding such 
a quantity for small systems is relatively routine on a classical computer. 

Feynman in fact analyzed two alternatives for simulating physics with computers ||. One 
uses a probabilistic classical computer that would produce from the same input as given to 
a physical system the same distribution of outputs as observed for the physical system. The 
other uses a computer constructed of distinctively quantum mechanical elements that obey 
the laws of quantum mechanics. This latter proposal is the quantum computer. 
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To the question, "Can quantum systems be probabilistically simulated by a classical com- 
puter?" , Feynman's answer was unequivocally "No." [] This answer is surprising for even at 
that time some quantum systems were being very successfully simulated probabilistically on 
classical computers, mainly by quantum Monte Carlo (QMC) methods ||. To the question, 
"Can quantum systems be simulated with a quantum computer?" , his answer was a qualified 
"Yes." He believed almost certainly that this could be done for a system of bosons but was 
unsure that it could be done for a system of fermions. In this paper we present a design for 
a universal quantum computer that will simulate a system of fermions. Before doing so, we 
first discuss some problems that can be solved by a probabilistic simulation of a quantum 
system on a classical computer and others that cannot. 

Probabilistic simulations of quantum systems on a classical computer are mainly per- 
formed with the use of the Monte Carlo method. These statistical approaches were in- 
troduced to overcome the difficulty of exponentially growing phase spaces by numerically 
evaluating the accompanying many-dimensional integrals by sampling from a function as- 
sumed to be non-negative. On a classical computer one can probabilistically simulate a 
quantum system like liquid He 4 |5| and produce results that accurately compare with ex- 
periment. The situation, however, is far from satisfactory. An unsatisfactory state of affairs 
results from the frequent breakdown of the non-negativity assumption and is called "the sign 
problem." The sign problem is manifested by the seemingly exponentially hard task of esti- 
mating the expectation value of an observable with a given error. Interestingly, Feynman's 
negativism about quantum systems being probabilistically simulated by classical computers 
was a claim that negative probabilities were unavoidable because of the "hidden variable" 
problem and the possible violation of Bell inequalities. The extent to which the sign prob- 
lem is a hidden variable problem is unclear. On the other hand, QMC methods do not 
faithfully adhere to Feynman's idea of a probabilistic computer. Two important differences 



There is as yet no mathematical proof that this is the correct answer. 
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are that most QMC simulations are non-local and performed in imaginary time. Feynman 
discussed real-time simulations on a local computer. Implications of these differences have 
been noted by Ceperley ]7) who suggests Feynman really argues just against simulating 
quantum dynamics on a local classical computer. In any case, probabilistic simulations on 
a classical computer clearly do not qualify as a universal computational scheme for general 
quantum many-body problems. The limiting factors, for whatever reasons, are negative or 
complex-valued probabilities whether the simulations are done in real or imaginary time. 

To place the sign problem in a better perspective, we will start with a real-time analysis of 
a collection of interacting quantum particles. Quantum mechanics tells us that these particles 
either obey Bosonic statistics, whereby the wave function is symmetric with respect to the 
exchange of the states of any two particles, or obey Fermionic statistics, whereby the wave 
function is antisymmetric (changes sign) with respect to the exchange of any two particles 
||. Examples of bosons are photons and gluons; examples of fermions are electrons, protons, 
neutrons, and quarks. Often these two quantum statistics conveniently and efficiently map 
onto a third, quantum spin statistics ||. Still in other cases, when particle exchange is 
unlikely, particle statistics is simply ignored. 

For a given initial quantum state |\l/(0)), a quantum computer solves the time- dependent 



Schrodinger equation 



ih^ = H\V) (1) 



by incrementally propagating the initial state via 

\%(t)) = e-* AtH / h e - iAtm ■ ■ ■ e- iAtH/h y \ty(0)) . (2) 

M factors 

(t = MAt and the Hamiltonian H is assumed time independent). It should be reasonably 
apparent that if the Monte Carlo method is applied to the evaluation of the right-hand 
side of this equation, it is faced with sampling from oscillatory integrands that are not 
always positive and have unknown nodal surfaces. Further, as time t increases, the integrand 
fluctuates with increasing rapidity. While clever stationary-phase forms of the QMC method 



have been developed, acceptable solutions are possible only for relatively short times. This 
form of the sign problem is called the dynamical sign problem, and we are unaware of any 



efficient |TIJ real-time QMC simulations for bosonic, fermionic, or quantum spin systems. 

Years ago ||, before quantum computers were proposed, it was realized that by trans- 
forming Schrodinger's equation to imaginary-time r via t — > — ihr the problem with the 
rapid fluctuations was eliminated. With this transformation, called Wick's rotation, one 
solves the diffusion-like equation 



dr 

by incrementally propagating the initial state via 



H\V) (3) 



gr( r )) = f - AT V ATjf • . • e - ATH j¥(0)) . (4) 



"V" 

M factors 



(r = MAt and the Hamiltonian H is assumed time independent.) This transformation 
permits QMC simulations of time-reversal invariant interacting boson systems to a high 
degree of accuracy. For systems of interacting quantum spins and fermions (or bosons with 
complex hermitian Hamiltonians |TI|]), the transformation does not solve the sign problem. 
For quantum spin systems, the difficulty is finding a basis in which all matrix elements of 
the positive-definite operator exp(— AtH) are positive. Most notably this difficulty occurs 
for frustrated quantum spins. For fermion systems, the problem is the Monte Carlo process 
causing state exchanges that because of the anti-symmetrization requirement just happen 
to produce samples which are as frequently positive as negative. For the sign problem found 
in both types of systems, the statistical error of the measured observables grows exponen- 
tially fast with increasing system size. Another difficulty with the imaginary-time approach 
is analytically continuing the results back to real-time if real-time, i.e., truly dynamical, 
information is needed [|I!J. This continuation is an ill-posed problem whose solution places 
extraordinary demands on the simulation ]TJ| . 

In this paper, we will focus on the dynamical sign problem for a system of fermions, 
seemingly the most challenging case. Eventually we will give a detailed implementation of a 



simulation of the dynamical properties of a collection of interacting fermions on a quantum 
computer. The implementation avoids the sign problem. First, in Section [TI| we will discuss 
more fully the mathematical origin of the dynamical sign problem in classical computation 
and show why a quantum algorithm overcomes the problem. In Section [TIT] we will give 



the elements required for Deutsch's quantum network model of a quantum computer flU 



The quantum gate in this model conveniently allows the propagation of systems of local 
two state objects, e.g., a localized quantum spin-| particle called qubit. We also propose 
a universal set of quantum gates (unitary operators) that allows generic propagation of 
systems of fermions (the fabled "Grassmann chip" ||15|| ). The resulting fermion algebra has 
been the main technical tool for studying the classical Ising model in two spatial dimensions 
16| , a prototype lattice model that had an enormous impact on our understanding of phase 



transitions. Next, In Section IV, we show how this propagation can be effected by the 
quantum spin gate. We will demonstrate the polynomial scaling of the construction of the 
initial state, its subsequent time propagation, and the measurement of some observable. 
Here we will also demonstrate the control of the error in the results. In Section |V], we apply 
our model of dynamical fermion computation to a toy problem to illustrate our procedures in 
more detail. The exact solution of this toy problem will be given in the Appendix; however, 
if one were to solve this problem in the obvious Hilbert space, its solution would require a 



quantum computer for a sufficiently large system. Finally, in Section |\T], we summarize and 
make some remarks about future research directions. 

Our universal fermion gate and its mapping to the standard universal quantum gate is 



very similar to the one recently discussed by Bravyi and Kitaev |L7]] who actually propose 
that a quantum computer built from fermions might be more efficient than one built from 
distinguishable two state systems. 
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II. DYNAMICAL SIGN PROBLEM 



In order to understand the mathematical origin of the dynamical sign problem we use 



the Feynman path integral formulation for continuum systems in the first quantization 
representation. In this formalism one maps a quantum problem in D dimensions into a 
classical one in D + 1 dimensions and then simulates that problem probabilistically on a 
classical computer. The algorithm is efficient except for the repetition needed to obtain 
sufficiently good statistics. The "distinguishable particle" quantum mechanical propagator 



of a system represented by the Hamiltonian H — Yli=iPi + ^(^) * s expressed as |L9 



G{Ti -> K'] t) = (W, t\e- iHt \n, 0) = / V[K{t)] e iS[n{t)] , (5) 

where the measure T>[TZ(t)} = YuxiM^oc^'Kit/ M)~ MD I 2 dlZ\ ■ ■ ■ (ITZm-i, and the action 



r-n(t)=-R.' 



1 fdK(r) 



o 



S [R(t)} = I dr{-\ ) - V(TZ(r)) } . (6) 



Bosonic or fermionic statistics are introduced by applying the corresponding symmetrization 
operator to the propagator, Eq. |^. However, because the dynamical sign problem occurs for 
any particle statistics, we will ignore particle statistics for the sake of simplicity. 

The description of the properties of different physical systems in terms of correlations 
of physical observables is the natural way to compare with available experimental informa- 
tion. In this regard, linear response theory provides a way to compute the response of a 
system to a weak external dynamical perturbation f20fl . This linear response is always ex- 
pressed in terms of a time correlation function of the dynamical variables that couple to the 
perturbation. For example, if we were to apply an external time- dependent magnetic field 
and we wanted to calculate the average induced magnetization, we would have to compute a 
time-dependent magnetization-magnetization correlation function. The two-time correlation 
function between arbitrary local dynamical variables A and B is 

C AB (t) = (A(t)B(0)) = (e iHt Ae- iHt B) , (7) 

if the Hamiltonian is time independent. Expressed in integral form 



_ / dKdK'dK" p(K, K")G(K -> K"; -t)A(K')G(K -> 72.'; t)5(72) 
AsW ~ / dKdK'dK" p{K, K")G{K' -»• 72"; -»• 72'; t) ' ^ 

where the "distinguishable particle" density matrix p specifies the way the system was 

initially prepared. The next step consists of writing Eq. |8| in path-integral form using Eq. |5| 

_ / dK x ---dK 2M AjKM^BjKi) P(K 1} - ■ ■ } K 2M ) e^-j^j (q) 
AB[) JdKi--- dK 2M P{K U ■ • • , K 2M ) e^-^M) ' > 

where P and $ are real- valued functions. Generically, a stochastic estimate of Cab(£) is 

° AsW " E { ^ } ~ (e»(W})) P ' UUj 

where the configurations {72j} are sampled from the probability distribution P (positive 
semidefinite measure). One immediately sees that the origin of the dynamical sign prob- 
lem is the oscillatory phase factor that leads to large phase fluctuations at long times. 
Manifestly, \{&®^ li ")p\ — >• in an exponential fashion as t gets larger. Therefore, the total 
statistical error for the evaluation of Cab(£) grows exponentially with time because of large 
cancellations both in the numerator and denominator. The so-called "fermion sign problem" 
is a particular case of this problem when e** = ±1 and time is imaginary [j21{| . 

Will a quantum computer solve this problem? At first glance it appears that it will. A 
quantum computer is a physical system, whether a system of fermions or not, and physical 
systems have no dynamical or fermion sign problems. Furthermore it has been argued that 
there are means for mapping most physical systems to a quantum computer in such a way 
that the quantum computer's controlled evolution mimics that of the physical system |3|,|22 



A closer look, however, makes the situation less clear. A quantum computer is a computer, 
and as such it suffers from limited accuracy. More importantly this type of computer pre- 
dicts results stochastically, meaning each measurement is one member of the ensemble of 
measurements possible from a distribution specified by the modulus squared of the wave 
function for the Hamiltonian H modeled by the quantum computer. For a fixed physical 
time t > 0, how accurate is an individual measurement, how accurate is the expectation 
value of these measurements, and how controlled is their estimated variance? Is the level of 
accuracy and control achievable polynomially with complexity and tl 
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There is an area where a problem similar to the sign problem has been recognized and 
resolved by quantum computation. Recently it was shown that quantum computation is 
polynomially equivalent to classical probabilistic computation with an oracle for estimating 
the value of simple sums of rational numbers called quadratically signed weight enumerators 
(QWGTs) |2"Bfl . In other words, if these sums could be evaluated, one could use them to 



generate the quantum statistics needed to simulate the desired quantum system. More 
specifically, what was demonstrated was the obtainability of expectation value of operators 
in quantum computation by evaluating sums of the form 

S{A,B,x,y)= (-l^s'V-' 61 , (11) 

b:Ab=0 

where A and B are 0-1-matrices with B of dimension nxn and A of dimension m x n. The 
variable b in the summand ranges over 0-1-column vectors of dimension n, h T denotes the 
transpose of b, \b\ is the weight of b (the number of ones in the vector b), and all calculations 
involving A, B and b are done modulo 2. The absolute value of S(A, B, x, y) is bounded by 
(\x\ + \y\) n - Quantum computation corresponds to the problem of determining the sign of 
S(A,\t(A), k,l) with the restrictions of having dg(A) = I, k and I being positive integers, 
and |5'(A,lt(A), k,l)\ > (k 2 + Z 2 ) n / 2 /2. dg(A) is a diagonal matrix formed from the diagonal 
elements of A and It (A) is a lower triagonal matrix formed from the lower triangular elements 



of A. Details of this quantum algorithm can be found in ||23|| . 

The main point is that these sums are similar to the numerator of Eq. [H| and attempts 
to estimate them by random sampling result in exponentially bad signal to noise ratios. 
In the case of QWGTs, quantum computers can estimate the sum exponentially better 
than classical computers, but the estimate is not exact. The situation for the dynamical 
sign problem is similar: Quantum computers cannot obtain exact values for the desired 
correlation functions, but can obtain estimates sufficiently exact to avoid the sign problem 
suffered by the known classical algorithms and to yield usable information about the physical 
models simulated. 

In this paper we will show explicitly how the sign problem is avoided in the case of 



simulating fermions. Below we will give a means for translating local fermion Hamiltonians 
into the Hamiltonians available in the standard model of quantum computation. In contrast 
to quantum simulations on a classical computer this translation prevents uncontrolled ex- 
change processes that are the dominant source of the fermion sign problem. With respect 
to the dynamical sign problem, we then argue by using standard error correction analysis 
developed for the standard model of quantum computing that these gates will enable suffi- 
ciently accurate measurements of correlation functions so the accuracy of the average of these 
measurements will be dominated by the statistical error. The statistical error is problem 
dependent but polynomially bounded, so that the difficulty associated with phase-weighted 
averages is eliminated. 

III. MODELS OF QUANTUM COMPUTATION 

The quantum control model of quantum computation assumes the existence of physical 
systems that can be controlled by modulating the parameters of the system's Hamiltonian 
Hp. The control possibilities are abstracted and used to implement specific quantum gates 
that represent the unitary evolution of the physical system over a time step obtained by 
specific modulations of the Hamiltonian. In most treatments, the physical systems, together 
with the gates, are then taken as the abstract model of quantum computation. The quantum 
control and quantum gate viewpoints are effectively equivalent, but to tie the computational 
model to the physics simulation problem more closely, we choose to describe quantum com- 
putation from the point of view of quantum control; that is, we will assume an Hp. In this 
context we begin by giving the standard model of quantum computation and then defining 
an alternative model based on fermions. 

Defining a model of quantum computation consists of giving an algebra of operators, a 
set of controllable Hamiltonians (Hermitian operators in the algebra), a set of measurable 
observables, and an initial state of the physical system. In the simplest case, the observables 
are measured by the method of projective measurements, and the initial state of the physical 
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system is an expectation value of the algebra's operators. 



A. Standard Model of Quantum Computation 

The standard model of quantum computation (Deutsch's quantum network representa- 
tion) is based on an assembly of two state systems called qubits, controlled by one- and 
two-qubit Hamiltonians, and on a measurement process determined by one-qubit observ- 
ables. 

Operator algebra: It is convenient to define the standard model through the algebra of 
operators acting on the qubits. This algebra is generated by the unit and Pauli matrices a x , 
a y and a z for each qubit j, 

i= C °) ; ^C o) ; ^0 o ! ) ; ^C -J- (i2> 

These matrices represent quantum operators with mixed commutation relations and span 
the space of complex-valued 2x2 matrices. For qubits j ^ k, the ex's commute, and for 
qubits j = k, they satisfy the relation a fJi iJ u + ovct m = 28^1, (/i, v = x, y, z). For a quantum 
register with n qubits, one may take the operator defined in terms of a Kronecker product 

o£ = 1 ® 1 ® • • • ® o> ® • • • ® 1 

j th factor 

of matrices acting on n two-dimensional linear spaces. Thus admits a matrix represen- 
tation of dimension 2" x 2 n . 

Control Hamiltonians: Control of qubits is achieved by applying Hamiltonians that act 
on either one or two qubits. A theorem P^LES! in quantum information processing says that 
a generic operation on a single qubit and any interaction between two qubits is sufficient for 
building any unitary operation. We take 

3 hJ 

where the a M (t) are controllable. Ideally, no constraints on the control functions are assumed. 
However, it is often simpler to design the required control by assuming that only one of the 
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ctn(t) is non-zero at any time. A quantum algorithm for this model of quantum computation 
consists of prescribing the control functions ||26|| . A convenient measure of the complexity 



of such an algorithm is the integral J Q dt' y^2 IM o^ t {t / ) (the action of the algorithm). The 
quantum gates are simply specific unitary evolutions that may be implemented in terms of 
Hp. A convenient universal set of gates is given by operators of the form exp(z<r^7r/4) and 
exp(icr*(j^7r/8). In the quantum network representation of the standard model, an algorithm 
is a specific sequence of these operators applied to the initial state of the qubits. 
Initial state: The initial state of the qubits is assumed to be an n term Kronecker product 
of the state |0) = Q which is an eigenstate of a 2 with eigenvalue 1. The state is completely 
determined by the expectation values (0|cr^.|0), which are 1 if the are all a l z or the 
identity, and are otherwise. Physically, the initial state has all "spins" up. 
Measurement: The final feature of the model of computation is the specific means for 
extracting information after a sequence of operations has been applied to the initial state. In 
the standard model, it is always possible to apply a projective (von Neumann) measurement 
f2?|j using the observables a\. With this capability, it is unnecessary to give an initial state 
explicitly, as the desired state can be prepared by using measurement and operations. To 
learn the expectation of an observable at the end of an algorithm, one repeats the algorithm 
and measurement procedure many times and averages over the measurements until the 
desired accuracy is achieved. 

For a description of the standard model of quantum computation in terms of quantum 



Turing machines, see ||28|| . Quantum networks are discussed in [24|. Introductory descrip- 



tions of the standard model may be found in [29,30 



B. Fermion Model of Quantum Computation: Grassmann chip 

Somewhat analogously, we now describe a standard model of fermion computation. For 
simplicity we only consider spinless fermions, i.e., fermions without internal spin degrees of 
freedom, although we could have considered more general fermionic algebras with internal 
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degrees of freedom ||. Physically, a system of spinless fermion might be a system of spin-| 
electrons in a magnetic field sufficiently strong to polarize it fully. The basic system of this 
model is a state (or fermionic mode) that can be occupied by or 1 spinless fermion. We 
define the model for n such modes. 

Operator algebra: We define the model through the algebra of the spinless fermion op- 
erators a - and at for each qubit j (J = 1, • • • ,n), i.e., through the algebra of 2n elements 
satisfying canonical anticommutation relations 

{a^cij} = , {a^a)} = 5^ , 

where {A, B} = AB + BA denotes the anticommutator or Jordan product, at (a •) creates 
(annihilates) a spinless fermion in state (mode) j. Each element admits a matrix representa- 
tion of dimension 2™ x 2™. The fermion algebra is isomorphic (as a *-algebra) to the standard 
model (or Pauli) algebra. The isomorphism is established through the Jordan- Wigner map- 
ping |g. 

Control Hamiltonians: We take 

Hp = ^2 \ a j(' t ) a j + <%(^) a j] + y^uijit) ( a i a i + a ) a i) ■ 

j ij 

This is a universal Hamiltonian, i.e., any other Hamiltonian for a system of interacting 
spinless fermions can be generated by it. Physical operators must be (Hermitian) products 
of bilinear combinations of the creation and annihilation operators, i.e., products of sums of 
ata 4 for arbitrary i and j modes. 

Initial state: The initial state is assumed to be an n term Kronecker product of the state 
|0) which is an eigenstate of the number operator ata^ with eigenvalue 0. The state is 
completely determined by the expectation values (0|ata-|0) = for all j. Physically, the 
initial state has all modes unoccupied. 

Measurement: Measurements can again be performed by using von Neumann's scheme of 
projective measurements. In Section [IV (J| , we will discuss another scheme more appropriate 
for the physical systems at hand. 
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In the next subsection we show how to simulate the fermion model by using the standard 
spin-| model. In particular it is possible to efficiently map the fermion Hamiltonians to 
Pauli operators which can be simulated using the control Hamiltonians of the standard 
model. This establishes that these two models of computation are polynomially equivalent. 
Here the point of view is similar to the one used for classical models of computation: the 
simulation of one model by another establishes their equivalence. 



IV. FERMION COMPUTATION VIA THE STANDARD MODEL 



In the previous Section we gave the elements required for Deutsch's quantum network 



model of a quantum computer |L4[] and proposed a universal set of quantum gates (unitary 
operators) that allows generic propagation of systems of fermions (the fabled "Grassmann 
chip" ||15|| ). Here we show how this propagation can be effected by the quantum spin gate. We 
will demonstrate the polynomial scaling of the construction of the initial state, its subsequent 
time propagation, and the measurement of some observable. We will also demonstrate the 
control of the error in the results. 

The first step is the observation that the set of 2n matrices 7^ (of dimension 2 n x 2 n ) 
satisfying the Clifford algebra identities 

{in, lv] = 25^ 

admits a representation in terms of Pauli matrices (Brauer-Weyl construction) 

7i = <?l, 72 = o\ 
7s = 74 = o\o\ 

7s = o*o%al, 76 = o\o 2 z ol 

n— 1 n—1 



i=i 

The following mapping of fermion operators 



3 ]a n 

3=1 i=i 
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0-1 



a i - ( 11-^ ) = ^•••^rvi = (-lr 1 72 ^ * 72j 



vt=l 



2 



a; -> | || -a z j (7+ = ^-ij- • • • a. yr 

vi=l 



where a± = x ™ v , defines a *-algebra isomorphism to the algebra of operators of the 
standard model. It is the so-called spin-| Jordan- Wigner transformation |JT|, and has the 
property that fij = a^a^ = a 3 + a J L. We note that fij is a "local" particle number (or density) 
operator and many types of interaction in physical systems are of the form "density times 
density" which simplifies the simulation as we will see. 

It is important to emphasize that the success of our approach depends upon the mapping 
of algebras (and not of Hilbert spaces). In this regard it is relevant to mention that the 
transformation just presented is a particular case of a more general set of mappings that we 
would like to name generalized Jordan- Wigner transformations |J . It is possible to imagine 
a quantum computer implemented, for example, with a three state unit (5*=1) instead of a 
qubit. In such a case, these generalized transformations still allow one to simulate fermions 
or particles with arbitrary statistics. 

Two additional comments are in order: The mapping for and a] described above 
corresponds to a one- dimensional array of spins. The extension to higher spatial dimensions 
can be done [32|53|| in various ways. A straightforward extension to two dimensions is to 
re-map the sites of a two dimensional array onto a one dimensional string and proceed as 
before. Also there is nothing special about using the fermion instead of a quantum spin 
as an alternative model of computation. One could have just as well used the hard core 
boson H . The main question is whether different algebras admit a physical realization. For 
hard-core bosons this realization is He 4 atoms. 
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A. Evolution 



Given a fermion model algorithm, it is necessary to efficiently obtain a corresponding 
standard model algorithm that at least approximates the desired evolution. The general 
principle is to map the time-dependent fermion Hamiltonian H{t) = J2i Hi to the standard 
model operators via the Jordan- Wigner transformation, express the result in terms of a sum 
of simple products of Pauli operators, and then use the Trotter approximation 



Each time step At is chosen so that the final error of the simulation is sufficiently small. 
Provided that the number of terms in the sum is polynomially bounded in the number n of 
qubits or fermionic modes and provided that each term can be polynomially simulated, the 
simulation is efficient in n and 1/error. 

To see how to do the simulation, consider the example of the bilinear operator H c = 
a-,a\ + clm\ in the control Hamiltonian of the fermion model: 



It is readily checked that the Jordan- Wigner transformation for the other terms in the control 
Hamiltonians are also decomposable into sums of a few products of Pauli operators. 

The whole idea of a quantum computer is simulating the operations we want by using 
unitary matrices U = exp(—iAtH P /h). These unitary matrices, representing quantum gates, 
perform reversible computation and are case-dependent. For our particular case, we know 
how to simulate H = a\ in the spin-| case (it is directly implemented in the standard model), 
so we ask what set of unitary operations produce the evolution U = exp(—iAtH c /h). In 
other words, how do we write a U = U 1 . . .U k such that H c = WHU? Consider for example 
the Hamiltonian H x = ala 2 ■ ■ ■ The procedure is as follows: The unitary operator 



-iAt(H +H 1 +...)/h TT -i. 



AtH */ h + 0((At) 2 ). 



(13) 



H c = (-iy \a\a\ ■ ■ ■ or V + + a\a\o\ ■ • • a{~ V_] 
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takes a\ — > cr*, i.e., U\a\U\ = o\. The operator 

1 + io\o\ 
takes o\ — > o^o 2 . The next step is 

U 3 = e^ a * a * 

to take cry(j 2 z — » — o^ofcrf. By successively similar steps we easily build the required string 
of operators: o\a 2 z • • • <ji~ x <j\.. 
If j is odd, 

will take a^a 2 ■ ■ ■ a{^ 1 — > (— l^^lo^cr^ ■ • • a J z , where [y] is the integer part of m/l. The 
final operator 

U j+1 = e l ^y 

will bring the control operator to the desired one (up to a global phase (— l^ 1 ^ - ]): 

If j is even, we need an additional unitary operator that flips the first qubit's <r* into a cr*. 
This flip is achieved with the operator 

/e" l f \ 
U j+2 = = (8) 1 (8) ••• (8) 1 

V e l l ) 

that takes <7* — > cr*. 

Hence, to construct this non-local fermion operator from the standard model requires 
additional steps that are proportional to j . This number scales polynomially with the 
complexity so the construction is efficient if the standard model is efficient. 

The one and two-body nature of naturally occurring interactions means that a term in 
a second-quantized representation of a Hamiltonian only has one of two forms: either a\a^ 
or aja^a^. We just demonstrated how to handle the first case. The second case merely 
requires applying that algorithm twice. This squares the complexity. 
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B. State preparation 



In this Section we discuss the preparation of states of physical relevance. Clearly, the 
preparation of the initial state is a very important step since the study and efficiency of the 
given physical process one wants to simulate depends upon it. 

Consider a system of N e fermions and n operators at (single particle states). A generic 
iV e -particle state of a Hilbert space Ti,N e of antisymmetrized wave functions can always be 
expanded in terms of the antisymmetric states 

N e 

\<5>a)=i[ b l I™ ) ' 

where bj creates a state j and |vac) = |0) <8> |0) • • • (g> |0) is the vacuum state (i.e., 6-|0) = 0, 
Vj). The operator is in general a linear combination of a^'s, i.e., fej = J^™ =1 ajPy where 
Pij is some matrix and N e < n. 

fn\ 

The states |$ a ) (a = 1, • • • , ) in general form an overcomplete set of non- 

\N e J 

orthogonal states that span the whole Hn e , i-e., redundantly generate Hn e - They are known 



as Slater determinants pOfl . Typically, |$ a ) is the result of a self-consistent mean-field (or 
generalized Hartree-Fock) calculation. Even a Bardeen-Cooper-Schrieffer superconducting 
state, which does not preserve the number of particles, can be written in this way after an 
appropriate canonical transformation which redefines the vacuum state |33 . 



One can easily prepare the states |$ a ) noticing that the quantum gate, represented by 
the unitary operator 

when acting on the vacuum state, produces |0) up to a phase . Therefore, the successive 
application of similar unitary operators will generate the state up to a global phase. 

Except for very small systems the total Hilbert space is too large to be fully used (it has 
an exponential growth with increasing system size). In practice, one works in a subspace 
of Hn £ that closely represents the physical state one is trying to simulate. Generically, as 
initial state, we will consider a very general expression of a many-fermion state: 



N 
a=l 

where the integer N is a finite and small number. The state can be prepared efficiently (in 
N) by a number of procedures. We now describe one. 

To make the description simple, we will assume J2a=i l a <*l 2 = 1 an d i^al^p) — da/3, 
which is equivalent to requiring to be an orthonormal set and \^/{t = 0)) to be 

normalized to unity. With these assumptions the steps of the state preparation algorithm 
are: 

1. Adjoin iV auxiliary (ancilla) qubits, each in the state |0), to the vacuum of the physical 
system. The resulting state is 

|0) <g> |0) <g> • • • |0) ®|vac) = |0) a <g> |vac) (14) 

V v y 

N 

2. From this state generate Yla=i a a\ a ) ® l vac ) where \a) is an ancilla state with only the 
a qubit being |1). The procedure for generating this combination of states is described 
below. 

3. For each a — 1, . . . , N, conditional on the a qubit being |1), apply the state preparation 
procedure for |$ Q ). The resulting state is 

TV 

J]a Q |a) ® |$ Q ) (15) 

a=l 

4. From this state generate 

1 N 

—= ^2 a al°)a <S> \$ a ) + terms without |0) a (16) 

V N a =l 

This step will also be described below. 

The final step is accepted if a measurement shows all the ancillas being returned to |0) a . 
The probability of successful preparation is thus J2a=i l a a| 2 /^ = 1/N. Consequently, on 
average, iV trials will be needed before a successful state preparation. 
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The procedure to produce step 2 is most easily described by example. We will assume 
N = 2. The problem then is to generate a 1 1 10) ® |vac) + a 2 |01) ® |vac) from 1 00) <8> |vac). In 
what follows all operations will be only on the ancilla part of the initial state so we will not 
explicitly show the vacuum. We also note that one can always apply a rotation to a given 
qubit that will take |0) into x|0) +y\l) with \x\ 2 + \y\ 2 = 1. The steps of the procedure are: 

2.1 Adjoin an ancilla qubit |b) initially being |0). The initial state is now |0) <8> 1 00) . 

2.2 Conditional on |b) = |0), rotate the a = 1 qubit, and then conditional on the a = 1 
qubit being |1), flip |b): 

a;i|0)<g>|00)+j/i|l)<g>|10) (17) 

2.3 Conditional on |b) being |0), rotate the a = 2 qubit, and then conditional on the a = 2 
qubit being |1), flip |b): 

x lX2 \0) <g> |00) + x iy2 \l) ® |01) + Vl \l) ® |10) (18) 

2.4 Project out the states with |b) being |1): 

ariJte|01>+3/i[10> (19) 

The rotations are chosen so that a± = yi and a 2 = %iV2- 

For the explanation of step 4, we will display the physical states. The problem is: Starting 
with ai|10) <g> + a 2 |01) <g> |$ 2 ), produce ([16]) . 

4.1 Adjoin an ancilla qubit |b) initially being |0). The initial state is now ai|0) (8> 1 10) <g> 
|$i) +a 2 |0) ® |01) ® |$ 2 ). 

4.2 Conditional on |b) being |0), rotate the a = 1 qubit, and then conditional on the a = 1 
qubit being |1), flip |b): 

ai(xi|0) (8) |00) + y x |l> ® |10)) <g> + a 2 (ari|0) ® |01) + y x |l> ® |11)) ® |$ 2 ) (20) 
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4.3 Conditional on |b) being |0), rotate the a = 2 qubit, and then conditional on the a = 2 
qubit being |1), flip |b): 

aiXifcalO) <g) |00) + y 2 \l) <g> |01>) <g> \^) + a 2 xi(a: 2 |0) <g> |00) + y 2 |l) (g) |01>) <g> |$ 2 ) (21) 

4.4 Project out the states with |b) = |0): 

x 1 x 2 (a 1 \00) <g> |$!) + a 2 |00) <g> |$ 2 » (22) 

The rotations are chosen so that x\x 2 equals l/V^/V where N = 2. Comparing step 2 with 
step 4, one sees they are structurally identical, differing by the set of amplitudes generated 
and the complementarity of the subspaces selected for the final result. This latter difference 
in some sense makes one procedure the inverse of the other. For the case of N > 2, one 
simply replaces steps 2.2 and 2.3 and steps 4.2 and 4.3 by "do loops" over a from 1 to N. 

On average, the entire procedure needs iV trials before a successful state preparation. 
(In many cases, the other measurement outcomes can be used also to avoid too many trials.) 
Construction of the initial state thus scales as 0(N 2 nN e ) < 0(N 2 n 2 ) so unless the number 
of Slater determinants is exponentially large, general many-fermion states can be initialized 
efficiently. 

C. Measurement 

While there is a variety of physical observables one measures experimentally and calcu- 
lates theoretically, at this time it is difficult to demonstrate that they all can be computed 
efficiently on a quantum computer. Fortunately, we will now argue that one important 
class of observables, the temporal correlation functions Cab^), can be computed not only 
efficiently but also accurately. These functions describe the temporal evolution of some 
observable A(t) in response to some weak external stimulus that couples to the system's 
variable -B(O). They are at the heart of understanding, for example, the optical properties 
of materials. 

21 



The goal is to determine correlations of the form Cab(£) = (A(t)B(Q)) = {e lHt Ae~ lHt B) 
up to a sufficiently small statistical error. Clearly, measuring efficiently Cab is not pos- 
sible for an arbitrary A and B. One sufficient condition is that A and B are efficiently 
simulatable Hamiltonians. This observation is based on a method for determining Cab re- 
fined by Kitaev [35] and applied to the measurement of correlation functions by Terhal and 



DiVincenzo |36|] . Here we give a different method based on an idea given in |37] . 



A general principle that can be used to obtain Cab is to decompose the operator whose 
expectation needs to be determined , i.e., A(t)B(0), into a small sum of operators of a simpler 
form and measure each summand individually. Our method directly measures expectations 
of the form (U'V) when algorithms for implementing the unitary operators U and V are 
available. General correlation functions are then determined by decomposing operators using 
a unitary operator basis, for example the one consisting of products of Pauli matrices. 

The method for measuring (U'V) consists of the following steps: 

f. Adjoin via a direct product an ancilla (i.e., an auxiliary) qubit a in the state (|0) + 
| f ))/ y/2 with density matrix p a = (l+cr£)/2 to the initial state of the system described 
by the density matrix p. 

2. Apply the conditional evolutions U x = |0)(0| <g> U + |1)(1| ® 11 and U 2 = |1)(1| <g> V + 
|0) (0| ® 1 (U = UiU 2 ). The methods of |24|] may be used to implement these evolutions 



given algorithms for U and V. 

3. Measure 2cx+ = + ia^ = 2|0}(1|. This may be done by measuring and a* in 
sufficiently many independent trials of these steps. 

4. Given the initial density matrix p, the expectation 

K + itf) P »® P = 2Tr n+1 [f/t| )(l| f /p a «g> p ] (23) 
= TW|0>(1| ® C/Vpa® p] = TinpVp] = (tfV) p , (24) 

as desired. The statistical noise in the measurement of (WV) P is determined by that 
of two binary random variables and therefore depends only on the value of Tr n [C/^Vp], 
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which is inside the unit complex circle. As a result it is a simple matter to determine 
the number of measurement attempts required to achieve sufficient statistical accuracy. 

The procedure for measuring Cab(£) can now De summarized as follows: First express 

m A m B 

A = A(0) and B = B(0) as a sum of unitary operators A = ^ Aj and B = ^ By . 

i=i j"'=i 

A convenient unitary operator basis that works well for the local observables of interest 
consists of all the products of Pauli operators, as each such product is easily implemented 
as a quantum algorithm. Then, for each j and j' one uses the just described method with 
U = e iHt A)e- iHt and V = By to obtain (Aj(t)B f (0)). V may be implemented by simulating 
the evolution under H, applying By, and then undoing the evolution under H. 

An alternative approach to the measurement process is von Neumann's projection 
method. We sketch it here for completeness and comparison. In this approach we also 
add an auxiliary (ancilla) degree of freedom to the problem. Suppose that this extra qubit 
corresponds to an harmonic oscillator degree of freedom |e). Then, we consider the composite 
state 

|*)s® |e) , 

where \^)s = Ylj l<A?')s ^ s ^ e state of the system we want to probe and \e) t is the state 
of the harmonic oscillator in the coordinate (x-)representation. The corresponding state in 
the momentum (p-) represent at ion is denoted \e) t . 

Assume the observable (t-independent Hermitian operator) we want to measure is A. 
Then, we are interested in determining s^l^Al^s in an efficient way. Suppose that we 
know how to implement the unitary operation U s (t) = e~ lAt . Following Kitaev we want to 
implement the following conditional evolution 

t 

From the spectral theorem we can write A = J2j A? s(4>j\- Then, 

U \<j>j) s ® |0) = J2 U ® l e )* = J2 e ' iAjt > e )* = l^i)* ' 
t t 
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where |0)o is a state with (p = 0) zero momentum. Basically, the conditional evolution U is 
a momentum translation operator for the harmonic oscillator extra state. Finally, 



U\$) s ® |6) = A J ® 1^')* • 



Although the second measurement method is conceptually simpler, it requires approxi- 
mately implementing the ancillary harmonic oscillator, the conditional evolutions for many 
different choices of t, and a more complex analysis of the measurement statistics. The con- 
ditional evolutions can be simplified somewhat, and in special cases (such as as a subroutine 



of factoring) become very efficient — see [35| 



D. Measurement Noise Control 

The quantum physics simulation algorithm described above is approximate and the out- 
put is noisy. In order to properly use it, we need to have explicit estimates of the error e 
in the inferred expectations given the noise in the implementation. Furthermore, the effort 
required to make e small must scale polynomially with 1/e. There are three sources of error 
that need to be considered. The first is associated with intrinsic noise in the implementation 
of the gates due to imperfections and unwanted interactions. The second comes from the 
discretization of the evolution operator and the use of the Trotter decomposition. The third 
is due to the statistics in measuring the desired correlation function using the technique 
given above. 



1. Gate imperfections 

The problem of gate imperfections can be dealt with by using quantum error correc- 
tion |38|]39|] and fault tolerant quantum computation pp|- fH|| . According to the accuracy 
threshold theorem, provided the physical gates have sufficiently low error, it is possible to 
quantum compute arbitrarily accurately. The fault tolerant computation implements unitary 
operations and measurements on encoded qubits with overheads bounded by (9(log fc (l/e)) 
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for some k. This exponentially efficient convergence implies that the effects of physical noise 
can in principle be ignored. 

2. Discretization error 

A second type of error is the one introduced by the discretization of the evolution oper- 
ator. This discretization is very similar to the one used in classical simulation of dynamical 
quantum systems. It is possible to estimate the size of this error by a detailed analysis of 
the discretization. For example using the Trotter approximation 

e -i(H 1+ H 2 )At = e -i(ffi)At/2 e -i(/r 2 )At e -i(fTi)At/2 + £>((At) 3 ). 

The coefficient of (At) 3 ~ —i(H 1 + H 2 ) 3 /Q can be bounded by estimating the largest eigen- 
value of Hi and H 2 . 

3. Measurement statistics 

Our technique for measuring the correlation function (A(t)B(0)) requires measuring the 
expectations of unitary operators UjVj, associated with the implemented evolution. In most 
cases, the operators A and B are a sum of 0{mA,B) products of Pauli matrices, so that 
OirriA) C/j's and 0{m B ) V^'s are needed. This means that the expectation is a sum of 
OirnxmB) random variables 7\,y, where |r.,y| < 1. To assure that the statistical noise (given 
by the standard deviation) is less than e it suffices to measure each r^y O^rriAiTiB/e 2 ) times. 

V. EXAMPLE: RESONANT IMPURITY SCATTERING 

A. Formulation of the Physical Problem 

Our toy problem is a ring of n equally-spaced lattice sites on which spinless fermions hop 
to nearest neighbor sites or hop onto or from an "impurity" state. The length of the ring is 
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L = na, where a is the distance between sites. The system is described by the Hamiltonian 
(in second quantized form) 

n -rrn 

H = -Tj2(ck +1 + cl +l c l )+eb^b+-=J2^ b + b ^^ » ( 25 ) 
i=i * i=i 

As usual, 6's and c's are fermion (anticommuting) operators. The index i labels the lat- 
tice sites (Ri = ia is the lattice site position) and strict periodic boundary conditions are 
assumed, i.e., 

4+n = 4 • (26) 

We now imagine that the system is initially prepared in the zero temperature ground 
state of the ring in the absence of the impurity. Then, at time t — 0, a fermion is injected 
into the impurity state. After the system has evolved for some time t, we want to compute 
the probability amplitude that the evolved state is still in the initial state. The relevant 
quantity to compute is (h — 1 and t > 0) 

G(t) = (^(0)16(^(0)1^(0)), (27) 
b(t) = e mt b(0) e- mt , (28) 



where the initial state is the Fermi sea of iV e < n fermions 

iVe-l 

|*(0)> = |F5>= Ilil )- ( 29 ) 

i=0 

|0) is the vacuum of fermions and 

n 

c L = ^E e ^' c 5- (so) 

V j = i 

The wave number kj is determined from the periodic boundary conditions, c\ +n = q, which 
implies 

= — — ^ , with rij an integer . (31) 

Lj 

There is no unique way to choose the set of n/s. The common convention is to define the 
first Brillouin zone as 
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with k values uniformly distributed in this interval with spacing Ak = 2n/L. 



B. Quantum Algorithm 



We want to write a quantum algorithm that allows one to compute G(t) (see the Ap- 
pendix for the exact closed-form solution). To this end, we start by representing fermion 
operators in Eq. [A5| in terms of Pauli matrices. Because of the form of the hybridization 
term, Eq. [A4], a most convenient representation is the following 



a 



.1-2 



™n- 



■1) a z a z ---a z a_ 



t 1 2 

4 = -W 



-l)XV, 2 ---a>r 1 , 



(33) 
(34) 
(35) 
(36) 
(37) 



from which the following mapping results 



tfb = 2(1 + a l z ) 



2(1 + a 
1 
2 



,^1^2 i ^1^2\ 



Therefore, the Hamiltonian operator reads 



H 
~2 



71-1 

j=0 



n— 1 



V 



]l + e ^ + ^^^ 2 + ^(a^ + a^) 



8=0 



(38) 



An additional simplification can be introduced when one realizes that the structure of the 
observable to be measured is such that 



b(t) = e im b($)e 



iHtu r\\ „-iHt _ e iHt a l e ~ iHt 



(39) 



where H is given by 
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H = 2eal + 2£ ko a 2 z + ^{a l x a 2 x + affi) , 



(40) 



and, therefore, the "string" one has to simulate has length equal to two (it involves only 
qubits 1 and 2) 

A{t) = b(t)tf(0) = e tRt a l _e- iRt al . (41) 

If we were to transform H = UH P iW unitarily with U = YYj=i e lH p^ tj and n a finite integer 
(UW = 11) in such a way that both H P1 and H P2 are physical Hamiltonians, then the 
simulation would be straightforward. (We call this type of mapping a physical unitary 
mapping.) For our two qubit case, one can always perform a physical unitary mapping with 



U 



(42) 



H P1 = -{E- VA 2 + V 2 )a] + -{E + VA 2 + V 2 )a 2 z , 



(43) 



with E = 2(e + £ ko ), A = 2(e - £ ko ), and cos9 = l/Vl + 5 2 with 5 = (A + VA 2 + V 2 )/V. 

In general, such a constrained transformation is not easily realized and one performs a 
Trotter decomposition 



e iHt = 


■ e iHs 


M 


e iH z s e iH xy s + Q ^ 


M 


(44) 


where H — H z + H xy with H xy 


_ V 
2 


+ fyl) and time slice s = h- 


On the other hand, 



one can easily perform a physical unitary mapping for e lHxyS 

e iH xy s _ jj e iH P1 sjj^ 

where H P1 = \{<y x — a 2 ) and 

{J = e i Ti < >l e - i l< r l e - i l< T Wz . 

Finally, the "string" one has to simulate with the quantum computer is 

A(t)c[S(s)] M ai [S\s)] M a\ 
S(s) = e iRzS Ue iHpiS W . 

and G(t) = (A(t)). 
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(45) 



(46) 



(47) 



VI. CONCLUDING REMARKS 



We investigated the implementation of algorithms for the simulation of fermionic quan- 
tum systems, and gave an explicit mapping that relates the usual qubit of a quantum 
computer to the fermionic modes that we want to simulate. Our attention focused on the 
so-called sign problem. It is a problem appearing in attempts to simulate classically the 
dynamics of quantum systems. We reviewed the origin of this problem and showed how 
this problem is avoided in quantum computing simulation. The evolution of quantum com- 
puters are intrinsically quantum mechanical and this is the main difference with a classical 
computer that allows one to solve the sign problem. We studied sources of errors in a quan- 
tum computer, such as gate imperfections and the expansion of the evolution operator, and 
argued that they would not open a back door to a problem similar to the sign problem. 

We gave a very general definition of what a model of quantum computation is. In 
particular and because of our particular interest, i.e., the simulation of fermion systems, 
we described the standard and the fermionic models ("Grassmann Chip"). These are, of 
course, not the only ones. Isomorphisms of *-algebras allow one to introduce more "esoteric" 
models 0. Indeed, there is nothing special about the spinless fermionic model of quantum 
computation. One could have used a "hard-core boson" model which admits, in principle, 
a realization in terms of He 4 atoms. The key point is the implementation of the physical 
gates. 

Our effort focused on the simulation of the dynamics of fermionic quantum systems. 
However other problems can be of interest: the thermodynamic or ground state properties 
of a Hamiltonian. Even if one had a quantum computer, it is not clear how to use it to 
efficiently compute these quantities. On the other hand, at present, no proof exists showing 
that this is not possible. 

An approach that in principle could be used to compute the spectrum of a Hamiltonian 
H (e.g., the ground state) or expectation values of arbitrary observables is the adiabatic 
"switching on" in conjunction with the Gell-Mann-Low theorem j|5| of quantum field theory. 
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The idea simply consists of introducing a fictitious Hamiltonian 

H £ (t) = H + f € (t) Hi , (48) 

where both H and H 1 are time independent operators (H = H Q + H x ) and the scalar 
function f e (t) is such that lim t ^ ±00 f e {t) = and lim^o f e {t) = 1, for an arbitrary adiabatic 
parameter e. In other words, H e (t = 0) = H and H e (t = ±oo) = H . H is typically 
an operator whose spectrum is known, e.g., an arbitrary bilinear operator representing a 
mean-field solution of H and whose eigenstates can be straightforwardly prepared (let's call 
it |$o))- The Gell-Mann-Low theorem asserts that 

lim u(o -oo)iy = 



if the state whose limit one is performing admits a series expansion in a coupling parameter 
characterizing the strength of Hi. This formal device generates the eigenstate adiabatically 
connected to |<3>o)- The theorem does not guarantee that if |$ ) is the ground state of H 
then |\& ) is the ground state of H. If the conditions of the theorem are satisfied then 
computation of the spectrum of H is straightforward. To our knowledge this approach has 
never been implemented in practice. 

The work presented here is only a first step in a program investigating the simulation 
of quantum systems using quantum computers. We have given a rather explicit algorithm 
for a simple problem and we will increase the complexity of the problems in the work to 
come. An interesting problem would be to provide algorithms to test for superconductivity 
in systems such as the Hubbard model. Such simulations using classical computers cannot 
unequivocally answer this important question because of the sign problem, but a quantum 
computer could. 
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APPENDIX: TOY PROBLEM: EXACT SOLUTION 

We can rewrite our original Hamiltonian, Eq. |25], in the wave number representation: 
1. Kinetic energy: 



T = -T 5>?c m + c\ +lCl ) = -^E el(k '~ k)R ' ( e "° + ^ > ( Al ) 



r 

n 

1=1 k,k 

Thus, 



T = Y l S kl c\c ki , (A2) 



i=l 



where Sk = — 2T cos ka . 
2. Hybridization energy: 

V i=l k,i 

Therefore, 



Thus, in the wave number representation, the total Hamiltonian reads 

n—1 n—1 

H = E e * C W + ttfb+Vj2 {ifi + b ] c k ) 5 ki0 (A5) 

i=0 i=0 

The quantization of k quantizes the single particle energy spectrum Sk- For example, for 
an n = 4 site ring, the allowed k values and energies Sk are shown in Fig. |I[ 
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FIG. 1. Single-particle energy spectrum for an n = 4 site ring 
Its ground state with n e = 3 is 



\FS) = 4o4i c iU°) , with k = , fa = - , and fc 2 



7T 



The Heisenberg equations of motion are: 



7T 

2 



6 f = i [if, 6 f ] 



ci. = z 



#,4 



where [A, B] = AB — BA. These commutators equal 

[H, 6+] = e fot + y c t 4 Q 

Clearly, one can distinguish two cases: 
1) k + : 



= e** 4(0) = e" 
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-i2tT cos fca f 



*■ 4(o) 



or in matrix notation 



2) k = (S = -2T): 

- i fet = e 6 t + y c t 

6t\ /e \/\ /fet 
cj ~ l [v S ) ' [cl 
The solution of these first-order differential equations are 

m )^ t (< v \ . (" m 

4(f)/ \v s>/ \4(o) 

From elementary matrix algebra 

e V\ ( % —y\ { e ttXl \ f x y 



U(t) = exp[i t 



V So) \y x J V e^ 2 / x. 



e V 

where Ai 2 are the eigenvalues of ( ) with corresponding eigenvectors 

V So. 



, and x + y = 1. After elementary algebraic manipulations: 



e + S , / / e - £ X 2 



-y 



and 



V 

x = 



v /( Al _ £ )2 + V 2 



V/(A 2 - e) 2 + V 2 ' 
Thus, the operator b(t) is expressible as 

b(t)=U* 1 (t) b(0)+U* 12 (t) c (0) , 

where 
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U^(t) = x 2 e~ itXl + y 2 e~ itX2 (A21) 
U* 2 (t) = xy (e- itXl - e~ itX2 ) , (A22) 

and detW(t) = e ft ( Al+Aa ) = e ft < e+6l) >. 
Finally, G(t) is simply given by 

G{t) =U* n {t) . (A23) 
\G{t)\ 2 = x A + y A + 2x 2 y 2 cos (Ai - A 2 )t and is plotted in Fig. |. 



|G(*)| 




o.o 



e = 0, V = 1/2 
e = J8, V = 1/2 
e = -8j F = 4 7 
e = 0, ^ = 4 



t 



FIG. 2. |G(i)| 2 for different values of the parameters (£q = —2). 



Note that G(t) is independent of the number of fermions present in the initial state 
|\l/(0)). This is why this toy problem reduces to a two qubit problem. 
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